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We study the intrinsic spin Hall effect (SHE) based on the orbitally degenerate periodic Anderson 
model, which is an effective model for heavy fermion systems. In the very low resistivity regime, the 
magnitude of the intrinsic spin Hall conductivity (SHC) is estimated as 2000 ~ 3000/ie~^Sl~^cm~^; 
It is about 10 times larger than that in Pt. Its sign is negative (positive) in Ce (Yb) compound 
systems with (/^'^) configuration. Interestingly, the obtained expression for the SHC depends 
only on the density of conduction electrons, but is independent of the strength of the c-f mixing 
potential and the mass-enhancement factor. The origin of the huge SHE is the spin-dependent Berry 
phase induced by the complex /-orbital wavefunction, which we call the "orbital Aharonov-Bohm 
effect". 

PACS numbers: 72.25.Ba, 72.25.-b, 75.47.-m 
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I. INTRODUCTION 

Spin Hall effect (SHE) is a phenomenon that an ap- 
plied electric field induces a spin current in a transverse 
direction. It has been attracting a great deal of interest 
as a method for creating and detecting spin current. Re- 
cently, the SHE in metallic systems are intensively stud- 
ied due to the interest for both the unsolved origin and 
the possibility of an application to spintronics device [1- 
9] 

Recent intensive studies of the SHE in transition met- 
als was initiated by the observation of the huge SHC in Pt 
[6, 9]. To elucidate the origin of the huge SHE in tran- 
sition metals, theoretical calculations of intrinsic SHE 
have been performed intensively [10-13]. The intrinsic 
SHE occurs in multiband metals with strong spin-orbit 
interaction (SOI) independently of impurities, which has 
a close relation to the intrinsic anomalous Hall effect 
(AHE) in ferromagnetic metals [14]. In ref. [13], the 
authors have revealed that huge SHEs are ubiquitous in 
multiorbital c?-electron systems by calculating SHEs in 
various 4c? and 5d transition metals. This study succeeds 
in explaining sophisticated and systematic experimental 
studies by Otani's group [9]. Therefore, it is strongly 
suggested that the intrinsic mechanism is dominant in 
transition metals. 

The large SHE in transition metals is induced by the 
phase factor of the d-orbital wavefunction in the presence 
of the atomic SOI, which we call the "orbital Aharonov- 
Bohm (AB) effect" [15]. The intrinsic SHC is predicted to 
be simply proportional to the spin-orbit polarization at 
the Fermi level {I ■s)fj_. According to the Hund's rule, the 
SHC should be positive (negative) in transition metals 
with more (less) than half-filling. Moreover, occurrence 
of large orbital Hall effect (OHE), which is a phenomenon 
that large d-orbital Hall current is induced by the electric 
field, is also predicted theoretically in many transition 
metals [15]. These fact suggests that a very large SHE 
and OHE may appear in /-electron systems compared 
to that in d-electron systems, since SHE and OHE are 
proportional to {I ■ s)^ and I, respectively 



In heavy fermion systems, very large AHE appears un- 
der the magnetic field [16-19]: In clean heavy fermion 
systems, anomalous Hall conductivity (AHC) CTgjj is inde- 
pendent of p sufficiently below the coherent temperature 
To, whereas dg^ oc above Tq, which indicates that the 
intrinsic contribution is dominant in such clean samples. 
In ref. [20] , they studied the AHE based on the orbitally 
degenerate periodic Anderson model (OD-PAM), which 
is an effective model for heavy fermion compounds. The 
obtained general expression has succeeded in explaining 
the huge AHC observed in heavy-fermion systems. Con- 
sidering the close relationship between SHE and AHE, 
one might expect that huge SHE can be realized in heavy 
fermion systems. 

In this paper, we study the intrinsic SHE based on the 
OD-PAM. It is found that the huge SHE in heavy fermion 
systems originates from the "orbital AB effect", which 
is given by the spin-dependent Berry phase induced by 
the complex /-orbital wavefunction. In the low resis- 
tive regime, the SHCs in Ce- and Yb-compound systems 
are predicted to be about 2000 ~ 3000/ie~^ • ^l~^cm~^ 
in magnitude, which are one order larger than that the 
value observed in Pt. The sign of the SHC is negative 
(positive) in Ce (Yb) compound systems with (f^^) 
configuration, since the SHC is proportional to the spin- 
orbit polarization {I ■ s)^ [15]. The obtained expression 
for the SHC does not depend on the strength of the c-f 
mixing potential nor the mass-enhancement factor. The 
SHC in /-electron systems will be measurable by using 
recently developed fabrication technique of high quality 
heavy fermion thin film [21]. 

Recently, present authors have studied the extrinsic 
SHE based on the orbitally degenerate single-impurity 
Anderson model (OD-SIAM) [22]. Using the Green func- 
tional method, we have derived both the skew scattering 
and side-jump terms analytically. It is found that the 
side-jump term derived in the OD-SIAM has a great sim- 
ilarity to the intrinsic term derived in the OD-PAM: The 
SHCs are simply proportional to {I ■ s)^ and their magni- 
tude are almost the same in both mechanisms. In section 
IV, we discuss the relationship between the intrinsic and 
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the side-jump mechanisms. 



II. 



MODEL AND HAMILTONIAN 



In the present paper, we study the intrinsic SHE and 
OHE for both Ce- and Yb-compound heavy fermion sys- 
tems based on the OD-PAM. In these systems, the num- 
ber of /-electron or hole is unity, and the total angu- 
lar momentum J is 5/2 or 7/2. In the presence of the 
strong atomic SOI, the J — 7 /2 level is about 3000 
K higher than the J = 5/2 level. Therefore, we con- 
sider only J = 5/2 (J = 7/2) state in Ce^+ {Yh^+) 
ion with 4/^ (4/^'^) configuration. We note that I ■ s = 
\ [J{J -f 1) - L{L + 1) - S{S + 1)] is given as follows: 



I ■ s 



-2 for J ^5/2, 



r s = - for J = 7/2. 

2 ' 



(1) 



Here, we introduce the following OD-PAM Hamilto- 
nian, which had been used to explain the large Van-Vleck 
magnetic susceptibility [23] and the small Kadowaki- 
Woods ratio [24] in heavy fermion systems with orbital 
degeneracy. 



H = 



J2 ^kcl^Cka + E^flj^jfkM' 
fccr kM 

{^kMafkM^ka + VkM crC^a IkAI^ 



Mka 



(2) 



where, cj,^ is the creation operator of a conduction elec- 
tron with spin cr — zLl. fj,j^j is the operator of a /- 
electron with total angular momentum J — 5/2 (7/2) 
and ^-component M (-J < M < J) for Ce^+ {Yh^+). 
Ek is the energy for c-electrons, is the localized /- 
level energy, and U is the Coulomb interaction for /- 
electrons. VkMa is the mixing potential between the c- 
and /-electrons, which is given by [20] 



VkM(j — 



2 J 



-Vi^VfYanLYriek,M, (3) 



where, is the Clebsh-Gordan (C-G) coefficient and 
Yi"^{6k, 4>k) is the spherical harmonic function. Here, the 
C-G coefficient for / = 3 is given by (for 1=3) 



a{{7/2~Mcj)/7Y'^5ra,M-a/2 for J = 5/2, 



= {(7/2 + Afa)/7} 



1/2. 



for J = 7/2. 

(4) 



Here, the fc-dependence of Vf is neglected due to the 
small radius of the /-orbital wave function. We also 
neglect the crystalline electric field splitting of i?-^-level 
since its effect on the intrinsic Hall effect would not be 



essential [25]. Hereafter, we put J/ = 0; the effect of 
Coulomb interaction on the SHC will be discussed in sec- 
tion IV. 

From the expression of the C-G coefficient in eq. (4), 
we see that conduction electrons with ^-spin mainly hy- 
bridize with M = -5/2 [M = 7/2) for J = 5/2(J '= 
7/2), which is consistent with the Hund's rule: That is 
the spin and orbital angular momentum are parallel (an- 
tiparallel) for J = 5/2 (J = 7/2). We wiU show that 
the sign of the SHC is explained by the spin-orbit polar- 
ization [15]. In the present study, we neglect the effect 
of crystalline electric field on /-orbitals, since it is small 
due to the small radius of the /-orbital wave function. 
Hereafter, we put h = I. 

In Fig. 1, we show the band structure of 
OD-PAM given in eq. (2). Here, rep- 
resents the hybridization bands given by E^ = 
1 

2 



{Sk + Ef)± V(efc --B/)2 +4|V^/|2 . In this study, we 

assume the metallic state, where the Fermi level fi lies in 
the c-f hybridization band. In this figure, kp the Fermi 
momentum and A = E^ — fi. 




FIG. 1: Band structure of the OD-PAM given in eq. (2). 
Here, is the hybridization band. 



'-''kM 



Gi. 



'-'kM' 



Vk, 



FIG. 2: The diagrammatic expression for the Green function 
in eq. (6) [20]. 

Here, the conduction and /-electron Green functions 
for OD-PAM in the absence of the magnetic field are 
given by as follows [20] : 



\VkMo 



, (5) 



(6) 
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We note that G'j,^^{uj) = [20]. The diagrammatic ex- 
pression for eq. (6) is given in Fig. 2. G^^is the /- 
electron Green function without hybridization given as 



0J + n-Ef 



(7) 



Now, we consider the quasiparticle damping rate r(w), 
which is mainly given by the imaginary part of the /- 
electron self-energy, Yik{<^) in heavy fermion systems. 
In the dynamical mean- field approximation (DMFA), 
the self-energy is composed local /-Green function, 
TjY^k'^kMM'iio) = g{uj)SMM', which is diagonal with 
respect to M and is dependent of M in the orbitally de- 
generate case [24]. Hero. is the number of fc-points. 
Therefore, in the present study, we assume that F is 
diagonal with respect to M, and is independent of the 
momentum. Moreover, since /-electrons arc degenerate 
in the present model, we assume that Fm is approxi- 
mately independent of M and can be approximated as 
Tmm' = l^MM', where 7 is a constant. In this study, 
we perform a calculation of the SHC using this constant 
7 approximation. Then, the retarded (advanced) Green 
functions are given by 



\Vf\ 



oj + H-Ef + (-)i7 



G' 



0/R(A) 



(w) = {oj + n - + \ 



(8) 



III. 



CALCULATIONS OF SHC AND OHC 



In this study, we calculate ctsh based on linear response 
theory. According to Streda [26], the SHC at T = in the 
absence of the current vertex correction (CVC) is given 
by CTSH = o-sH + ^SH' where 



-"SH 



TI 



2nN I ^ y 



Attn 



N ^ 



dwTr 



J; 



cj=0 



(9) 



CAR 



J,VG 



(Rf^ A) 



(10) 



Here, (Tgjj and (Tg^ represents the Fermi surface term and 
the Fermi sea term, respectively. 

In the present model, the charge current operator is 
given by = — et)fe^, where — e (e > 0) is the electron 
charge, and 



(11) 



Next, we explain the Sz-spin current operator j^. In the 



present model, Sz is given by 



Mfe, 



(12) 



M 



where Sm = J2ma 1 ['^mo-]^- It is straight forward to 
show that Sm = -f- (f-) for J = 5/2 (J = 7/2). Then, 
the spin current = |^fc^j •^zj /2 is given by 



(13) 



In a similar way, the total angular momentum current 
operator, = Jzj /2, is given by replacing Sm 

in eq. (13) with M. Then, the the orbital angular mo- 
mentum current operator, = Jf^ — J^, is expressed 
as 

= E{^(^-^M)^cUM+h.c}.(14) 

Then, the orbital Hall conductivity (OHC) ctoh = 
( J^) /Ey due to the OHE is given by ctoh = + '^oH' 
where ctqjj and ctq^ are respectively given by eqs. (9) 
and (10) by replacing with . 

Here, we study the velocity given by the c-/ mixing 
potential VJcMct [20]: 



dkx 



d 

+ {VkMa(X*M,a) OlM,a 



(15) 



Here, is the anomalous velocity given by fe-derivative 
of the phase factor aM,a = exp {i (M — ^) in VkMcr- 
Figure 3 is a schematic view of the anomalous velocity 
v"" oc Vk(t>k- Since oc ky and thus ^kV'^{dek/dky) 7^ 
0, the anomalous velocity gives rise to the large SHE 
and AHE in heavy fermion systems. On the other hand, 
oc kx gives a normal velocity. In eqs. (9) or (10), the 
terms which contain single give rise to the SHC. 



Calculation of the Fermi surface term 



Here, we calculate the SHC by neglecting CVC accord- 
ing to eqs. (9) and (10), using eqs. (11) and (13) 
and are composed of the conduction electron term 



4 



k Fermi 
y r surface 




FIG. 3: A schematic view of the anomalous velocity u". 



dsk/dk^ = d^£k and the hybridization term d^Vu- Fig. 
4 shows the terms for ctsh in which J^, Jy is composed 
of zero or one (9m«^fc- Fig. 4 (a) gives large SHC since 
d^Vk includes the anomalous velocity in eq. (15). We 
note that the terms in 4 (b) that are composed only of 
dx£k ■ dySk vanishes identically. Moreover, there exists 
the terms that are proportional to d^jVkduVk, as shown 
in Fig6. In Appendix B, we will show that these terms 
are much smaller than the contribution by Fig. 4 (a). 
Therefore, we here focus on the terms in Fig. 4 (a). 

In this subsection, we derive the analytical expression 
for the Fermi surface term, since the Fermi surface term 
dominates over the Fermi sea term, as discussed in pre- 
vious studies [10, 11, 13, 27]. The Fermi sea term will be 
derived in section IIIB. 




FIG. 4: The diagrammatic expressions for ctsh. (a) The di- 
agrammatic expressions for the dominant terms, (b) The di- 
agrammatic expressions of the terms composed only of fl^Efe, 
which vanishes identically. 



According to eqs. (9), (11) and (13), the Fermi surface 



term (jgjj for Fig. 4 (a) is given by 



2t:N 



kMa 
dVkMa dSk 



vc 



Gf(0)pG°^^(0)-K, 



(16) 



Here, we confine ourselves to the case J = 5/2 state 
corresponding to Ce "^^-ion. In section IV A, we will dis- 
cuss the case for J = 7/2 state. Then, by using the 
following relationships 



^M2|Wp = ^(l + 16sin20) 

Ma 

Y,a^\VkMa\''=2\Vf\\ 



Ma 



Y,M(T\VkMa\^ = |V}f (1 -4sin2 



(17) 
(18) 
(19) 



Mo 



1.2 I jU2 



1 sin 6* sin (/) dsk dsk . „ . , 

= ITT sm^sm 0,(20) 



sin^ 9 ' dky 
eq. (16) is transformed as follows: 



dk 



CTSH 



2ttN 



^ ,2 1 9£k 

k 



dk ifM - EkY + 7' 



\Gm\\ 



(21) 



where k = | fc | . 

Here, we analyze eq. (21) when 7 is small enough: In 
this case, 



where ik = £k 



Gr(0) = iM-efc 
^ (m- 



\Vf\' 



11 -Ef 
k + iVcY^ , 



and Fr 



\Vf\ 



(22) 



^■7. Since 



7/(0;^ +7^) — 7r(5(a;) for small 7, we obtain the following 
relationship: 



1 



(a* - Ek 



r, 



■Si^l-Ek) 



(23) 



Substituting above equation into eq. (21), we obtain 
the following relationship for small 7: 

e 52^1^,(^_^^), (24) 



fSH 



27r7V 21^ k dk 



Now, we approximate the conduction electron as free 
electron. Then, (t|jj for J = 5/2 is given by 



^SH 



2127r2 



e 26 
27ra 21 



Nfs, 



(25) 
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where a is the lattice spacing and Nps represents the 
number of large Fermi surface. The first line in eq. 
(25) means that the SHC depends only on the density 
of conduction electron ric = /cp/STT^, except for Nps. 
This result suggests that SHCs in Ce-compound heavy 
fermion systems take similar large negative; values. The; 
second line in eq. (25) is obtained by putting fcp = n/a. 
When a = 4A, then e/27ra w lOOOfe-^fl-^cm. If we 
assume that Nps = 2 ^ 3, we obtain ctsh = 2000 ~ 
3000/ie~^rj~^cm~^ for Ce-compound system. Interest- 
ingly, the expression obtained above is independent of 
the strength of the c-f mixing potential. 

Next, we discuss the Fermi surface term for the OHC. 
By replacing with in eq. (9), o-qh can be calcu- 
lated in the same way as SHC. The obtained result is 



7 7 



(26) 



Thus, ctqjj shows a large positive value in Ce-compounds. 

is satisfied in tran- 



In contrast, the relation ctqjj ^ \af 



SHI 



sition metals since the SOI is weak and {I ■ s)^ <C 1 [13]. 



Calculation of the Fermi sea terms 



In this section, we derive the analytical expression for 
the Fermi sea term Ugy , and show that the Fermi surface 
term (7) dominates the Fermi sea term (//). 



According to eqs. 10, 11, and 13, the Fermi surface term crg^ for Fig. 4 (a) is given by 



^11 — 



47rA^ ^ 

kMa 



du) 



3(7 



2 V 2 



+ Sm 



dkx 



dky [ 



du 



Gf(c.)-G{»^^-(i?^A) 



kMa 



dkr. dk, 



5e./aG-(.)^,«^^^_^.,^^^aGr(.) 



(27) 



Using the relations in eqs. (17) - (19), and performing 
the M, cr-summations in eq. (27), it is transformed as 



TI 



— e 



4wN 
Im 



I/: 



7 r ^^ k dk 



du 



[J-oo [{iJ-ek){cJ-Ef + i^)-\Vf\^^ 



(28) 



To perform the cj-integration in eq. (28), we rewrite the 
integrand in eq. (28) as follows: 

{lo + ij.- £fe) {uj + fi- + I-,) - |V}f 
= {u: + ti-E++i^+){uj + ii-E^+i^-),{2<d) 



where 



7 / Ek-Ef 
2 V ^ E+-E7 



(30) 



Then, the w-integration in eq. (28) can be performed 
analytically as follows: 
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1 + -2m-z(7++7') 

[E+ -E^-i{j+- 7-)]' {Ei-^l- *7+) {E^-l^- il- 



■In- 



k - M - ^7 



[K -E^-^ (7+ - 7-)] 



(31) 
(32) 



Wc analyze cqs. (31) and (32) when 7 is small: Since 
Im[ln(a; ± 17)] = =F7r0(a;), the imaginary part of eqs. (31) 
and (32) is approximated as 



Im{eq.(31)} 
Im{eq.(32)} 



{Et-E-^)' 
(m - E-^ 
{Et-E^f 



(33) 
(34) 



for Ce-compounds, where the Fermi level lies under E-^" , 
as shown in Fig. 1. Substituting above equations into 
eq. (28), a^l^ and (TgH is given by 



lib 



52, 



2ttN 



! 1 dek -^S (/i - E^ ) 



e ^52 ^ldek2Tre{ti-Ei^' 



27riV^ 7 ' ^' fc afc (i;+ 



(35) 



■• (36) 



We will explain in Appendix A how to perform the k- 
summations in eq. (36). In case of \Vf\'^/{Ef — /z) 1, 
final expressions for and are obtained as 



^SH — ^Sn- 



(37) 
(38) 



Here A;;^ = 
1 , W 



and a^^ = 



dek 



■. Considering the relation 



Ekp in 



{t^-Eff 

Fig. 1, it is straight forward to show that A^^ ~ 1 up 
to 0((A^/|y^|)^). In this case, we obtain the following 
relationships for small 7: 



„IIb Ila 



(39) 
(40) 



Therefore, two Fermi sea terms and almost 
cancel, and as a result, the Fermi surface term crgjj gives a 
dominant contribution to the SHC [10, 11, 13, 27]. Note 
that the same relations also hold for the OHC, and the 
total OHC is mainly given by the Fermi surface term. 



IV. 



DISCUSSIONS 



A. SHC and OHC in Yb-compound system 

Now, we discuss the SHC for J = 7/2, which corre- 
sponds to the case in Yb-compound systems. 

To perform M, cr-summations, we use the following re- 
lations for J = 7/2: 

$]M2|V^fcM.P = ^(l + 3Osin2 0), (41) 



Ma 



2Vf, 



(42) 



Ma 



^M(7|14Ma|' = V/(l + 3sin2 0). (43) 



Ma 



By using the above relationships, we can perform the 
calculation of crgjj by following section HI A. As a result, 
cTgjj for J = 7/2 takes a large positive value as 



15 fcp 
^14 2^^^^ 



e 15 
2^14 



Nfs- 



(44) 



The second line in eq. (44) is obtained by putting fcp = 
7r/a. This result suggests that SHCs in Yb-compound 
heavy fermion systems take similar large positive values. 
We can also calculate the Fermi sea for J = 7/2 by fol- 
lowing section IIIB. Then, we recognize the relationship 
in eqs. (39) and (40) for J = 7/2. 

In the same way, the OHC for J = 7/2 state is given 
by 



(45) 



Therefore, we note that the sign of SHC is negative 
for J = 5/2, it is positive for J = 7/2, whereas the 
OHC is positive for both cases. These facts are consistent 
with the results obtained in 4d- and 5d transition metals 
[13, 15]. In section IV B, we will show that the sign of 
SHC is equal to the sign of the spin-orbit polarization 
{l-s)^[15]. 
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conduction 
electron 




orbital AB phase factor = e'^^ 



FIG. 5: Effective Aharonov-Bohm phase in two-dimensional 
OD-PAM. 



B. Orbital Aharonov-Bohm Phase Factor 



In previous sections, we have discussed the SHE based 
on the OD-PAM using the Green function method. In 
this section, we give an intuitive explanation for the ori- 
gin of the huge SHE in heavy fermion systems. For 
this purpose, we consider the two orbital model with 
M = ±5/2, assuming the strong crystalline field. In 
the case of J = 5/2, the c-f mixing potential is given 

by VMa{f) cx { 76^3^3" (r)5Af,-5/2. + Yi^r)SM..5/2a} in 
the real space representation: If we drop the second term, 
it can be approximated as VMair) cx Y^^'^{f) cx q-^^'^'P^ ^ 
where 0,. = ta.n^^{y/x). 

In Fig. 5, two examples of the clockwise motion of 
the conduction electron along the nearest three sites 
[/m — >■ c — >■ c — > /a/] are shown. Here, 9i represents 
the angle between the incoming and outgoing electron. 
Therein, the electron acquires the phase factor e"^"'^ due 
to the angular dependence of the c-f mixing potential 
in real space, Vm(t(t^)- This phase factor can be inter- 
preted as the "orbital AB phase factor" at the /-site, 
which works as the effective magnetic flux {—3a9/2Tr)(j)o 
through the area of the triangle. Here, (f>o = 27Th/e is the 
flux quantum. On the other hand, Vm^ is approximately 
given by VMa{r) cx e^*"^"*" for J = 7/2. In this case, the 
effective magnetic flux per triangle is {3a9 /2Tr)(f>Q, which 
is opposite to that for J = 5/2. 

In summary, a conduction electron acquires the spin- 
dependent "orbital AB phase factor" , which originates 
from the spin-dependent c-f hybridization in the pres- 
ence of strong SOI. This is the origin of the huge SHE 
in heavy fermion systems. This consideration also ex- 
plains the sign difference of the SHC between Ce- and 
Yb-compounds. Thus, the origin of the SHE in heavy 
fermion systems is well understood based on the simpli- 
fied two-orbital model. 



C. The relationship between the intrinsic and 
side-jump terms 

So far, we have studied the OD-PAM with translational 
invariance, and found that huge intrinsic SHC emerges. 
Here, we consider the depletion of /-electron. The quasi- 
particle damping rate 7 increases in proportion to the 
depletion ratio x. In the case of x ^ 1, the intrinsic SHC 
is independent of a: if 7 is smaller than the band splitting 
[11, 13]. In addition to the intrinsic term, the depletion 
may induce the extrinsic terms, that is, skew scattering 
term a^"" and side-jump term a^^. 

In the dilute limit where I — x 1, intrinsic term does 
not exist. In this case, present authors had studied the 
extrinsic SHE based on the orbitally degenerate single- 
impurity Anderson model [22]. For fci? = 7r/a (a is a 
lattice spacing), the expressions for skew scattering and 
side-jump terms are obtained as 



skew 



1 



= —S2-{l-S 



2iTa 
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(46) 
(47) 



for both J = 5/2 {{I ■ s)^ = -2) and J = 7/2 {{I ■ s)^ ^ 
3/2). Here, 62 is a phase shift for d partial wave. From 
the above equation, we find that the extrinsic term is 
proportional to the spin-orbit polarization {I ■ s)^. In 
these two Anderson models, both intrinsic term and 
side-jump term originate from the anomalous veloc- 
ity that arises from the ^-derivative of the phase factor in 
the mixing potential. Here, we compare eqs. (25), (44), 
and (46). Very interestingly, the following relationship 
holds in a accuracy of ±7.2%: 



'SH 



(48) 



This fact indicates unexpected close relationship between 
the intrinsic term and the extrinsic side-jump term, and 
therefore it would be very difficult to distinguish these 
two mechanisms experimentally. This fact would be the 
reason why intrinsic (or side-jump) term are widely ob- 
served from single crystals to polycrystal or amorphous 
compounds. 



V. SUMMARY 

In this paper, we studied the intrinsic SHE and OHE 
based on the OD-PAM. We derived the analytical ex- 
pression for the intrinsic SHC and OHC based on the 
linear response theory. Both SHC and OHC are mainly 
given by the Fermi surface term (/). The obtained re- 
sults for Ce-compounds (J = 5/2) are given by eqs. (25) 
and (26), and those for Yb-compounds (J = 7/2) are 
given by eqs. (44) and (45). The SHCs for both com- 
pounds are approximately expressed by eq. (46). These 
results suggests that SHCs in Ce- (Yb-) compound heavy 
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fermion systems take similar large negative (positive) val- 
ues; 2000 ~ SOOO&^^fi^^cm"^ in magnitude. The mech- 
anism of the huge SHE and OHE in /-electron systems 
is the "orbital AB effect", which is given by the spin- 
dependent Berry phase induced by the complex /-orbital 
wavefunction. Therein, the SHC is proportional to the 
spin orbit polarization {I ■ s)^. The SHC in /-electron 
systems will be measurable by using recently developed 
fabrication technique of high quality heavy fermion thin 
film [21]. 

Here, we briefly comment on the effect of the Coulomb 
interaction U. In the present study, we have calculated 
the SHC with U = 0. In the PAM, the effect of the 
self-energy correction is represented by the renormaliza- 
tion of the mixing potential VkAia — ^ ^/zVkMa, where 
= 1 — ^S(e) — m* /m is the renormalization factor 
due to the self-energy [28, 29]. Since the SHC obtained in 
this study is independent of VfcMo- , the SHC will be inde- 
pendent of the mass-enhancement due to Coulomb inter- 
action. (In contrast, the AHE under the magnetic field is 
proportional to the magnetic susceptibility c< m* /m.) 
Next, we discuss the CVC due to Coulomb interaction. 
In ref. [20], it was proved that the CVC by U does not 
give rise to the skew scattering term, and thus its quan- 
titative effect on the SHE is expected to be small [20]. 
However, the CVC due to spin fluctuations might be sig- 
nificant in nearly quantum-critical-point [30]. This is an 
important future issue. 
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where X = - - When \Vf\'^/{Ef - fi) > 1, 

the first term in the bracket in eq. (Al) is approximated 
as « 1. As a result, a^f^ is given by 



Appendix B: Calculations of the term proportional 

to dxVkdyVk. 




5 (j + Sm) 



FIG. 6: The diagrammatic expression for the term propor- 
tional to dxVkMadyVkMa- 

In the main text, we have calculated the term propor- 
tional to d^Vdv£k, and explained that it gives a domi- 
nant contribution to the SHC. In this appendix, we de- 
rive the SHC given by dxVkdyVk, and show that it is 
very small and negligible. In this case, to perform the 
M, (T-summations, we use the following relations: 



Appendix A: Derivation of eq. (38) 

Here, we explain the way we performed the k- 
summations in eq. (36), and derive eq. (38). In per- 
forming the fc-summations analytically, we assumed that 
the density of state N{ijj) for conduction electron is con- 
stant: Y.k = ^(0) ! <^^k- Then, 



X 



X 



dei 



{K-Ek 



d£k 



-co [{£k-Eff+A\V}\^\ 

I { X ^ 



3/2 



A2 



(Al) 



E 



I fa M\dVkMadV, 



2 V 2 



7 



kMa 



1 



dk„ 



12 

i-cos'.^,, 
^ 1 fa M\ dVkM^ dV*^,^ 



M,AI' 



2 V 2 



dk,. 



■16 1 



1 f a M\ dVkMa H. dVkM'a' 



M,M' ,a,<j 

16 



2 V 2 7 
1 



dkx 



dk,, 



(Bl) 



Here, we first perform the calculation for the Fermi 
surface term. By using the above relationship shown in 
eqs. (Bl), the SHC given by (A)-(D) in Fig. 6 is given 
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by 

I(A-D) 



2^7- 



^^|y/pIm{G^«(0)G=^(0)} (B2) 



Substituting the above equations into eqs. (B2) and (B3), 
and performing the fc-summation, we obtain 



+ 



e 64 
2^21 



5^ -i|F/rReG/«(0)ImG/«(0)|G^(0)p. 

(B3) 



The diagrammatic expressions for eqs. (B2) and (B3) are 
respectively given by (A) and (B), and (C) and (D) in 
Fig. 6. The contributions from the diagrams (E) and (F) 
turn out to cancel out. 

When 7 is small, we obtain a following relationship: 



Im{G^^(0)G"^(0)} 



7T5{fl- Ek) 



ReG/^(0)ImG/^(0)|G^(0)| 



2 1 TTd{n - Ek) 



(B4) 



\Vf\' H-Ef 



(B5) 



-"SH 



e 10 
- — —aNps, 
zira 21 



(B6) 



where a is defined by — Skp \ = ae^f 



In a similar way, we calculate and o-gn • After 
performing the M, cr-summations using eq. (Bl), we ob- 
tain the following expression of the Fermi sea term for 
(A)-(D) in Fig. 6: 



IHA-D) _ 



_e_48 
4^21 

e 128 
2^"2r 



(B7) 
(B8) 



As explained in section IIIB, after performing to- 
integration, above expressions are rewritten as follows 



J 



for small 7: 



Ila(A-D) 
''SH 



\Vf\' 



E+ - E; 



2n 21 ^ {Ef - E^)iE+ - E^y 



(B9) 



IIb{A-D) _ 
''SH 



■ 487r 



y- 



\Vf\' 



+ ■ 



{Ef - E^){E+ - E^f {Ef - E^){E+ - E, 



(BIO) 



Performing the fc-summations in eq. (B9), and as a re- 
sult, we obtain the following expressions for <J^^: 



Ila(A-D) 
''SH 



2 16 ■ 

— p a 

7^ 21 



Nfs, (Bll) 



where p is defined by \E'^^ — ^| = jjekp. As recognized 
in Fig. 1, the relation a « /3 ~ ^ is satisfied since E^^ w 
Ekp is satisfied in the present model. Since the relation 
Afep = 1 -h 0((A/y/)^) holds well as discussed in ref. 
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IIIB, (TgH is given by 



(7, 



IIa{A-D) _ 



e 10 



(B12) 



^" 27ra21' 
To perform fe-summations in eq. (BIO), we use the fol- 

fc2 ~ ^(0)]^ / dsk. Then, 



lowing approximation: ^ 
<^SH is given by 



e 
2'Ka 

e 1 
'2^42 



1 

14 



1 

21 



A 



1 14. 



Npsoi. 



Nfs 
(B13) 



Finally, we obtain the final expressions for o"|jj,(Tg^", 
and a^f^ arc given by the siimmations of eqs. (25) and 
(B6), eqs. (37) and (Bll), and eqs. (38) and (B13), 
respectively. 



^I-tot 



//a— tot 



^//5-tot 



■5^(S + S"I''- 



Here in eqs. (B14) -(B16), the terms that is proportional 
to a is given in Fig. 6. In total, the SHC is given as 



26 



1 



^^- = -2^ Ul + i^^'^^^- 



CBlT) 



In eq. (B17), the factor 26/21 and q:/42 in the bracket 
come from the terms with d^Vkdi^Sk and the terms with 
3/iVfe9yVfe, respectively. Since a ^ 1/2, the terms pro- 
portional to d^VkdySk shown in Fig. 4 gives a dominant 
contribution. 
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